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Abstract 

Goals, Scope and Background. Eutrophication and hypoxia, 
which are already serious environmental issues in the Midwest¬ 
ern region of the United States and the Gulf of Mexico, could 
worsen with an increase emphasis on the use of corn and soy¬ 
beans for biofuels. Eutrophication impacts from agriculture are 
difficult to integrate into an LCA due to annual variability in 
the nutrient loads as a factor of climatic conditions. This vari¬ 
ability has not been included in many relevant energy or row 
crop LCAs. The objective of this research was to develop a rela¬ 
tively simple method to accurately quantify nutrient loadings 
from row crop production to surface water that reflects annual 
variations due to weather. A set of watersheds that comprise 
most of eastern Iowa was studied. Ample data describing corn- 
soybean agriculture in this region and nutrient loadings to the 
Mississippi River enabled the development, calibration and vali¬ 
dation of the model for this particular region. 

Methods. A framework for estimating lifecycle inventory data 
for variable nutrient loading from corn-soybean agriculture was 
developed. The approach uses 21 years of county-average data 
for agricultural and annual rainfall for 33 counties that approxi¬ 
mate three major watersheds in eastern Iowa. A linear equation 
describes the relationship between the fraction of the applied 
nutrients that leach into the surface water and the annual rain¬ 
fall. Model parameters were calibrated by minimizing the error 
in the difference between actual and modeled cumulative dis¬ 
charge to the Mississippi River over the period 1988-1998. Data 
from 1978-1987 were used to validate the method. Two sepa¬ 
rate approaches were then used to allocate the nutrient flows 
between the corn and soybeans. 

Results and Discussion. The total nitrogen (TN) and total phos¬ 
phorus (TP) leaching models provide good representation of the 
variability in measured nutrient loads discharged from eastern 
Iowa watersheds to the Mississippi River. The calibrated model 
estimates are within 1.1% of the actual 11-year cumulative TN 
load and 0.3% of the TP load. In contrast, a standard method 
used in other lifecycle assessments for estimating nutrient leach¬ 
ing based on a constant fraction of the nutrients leached pro¬ 
vides a reasonable average, but does not capture the annual vari¬ 
ability. Estimates of the TN load that can be allocated to corn 
range from 60 and 99% between two allocation methods con¬ 
sidered. This difference stems from a poorly understood sym¬ 
biosis of nitrogen flows within the corn-soybean rotation that is 
difficult to integrate into an LCA. 

Conclusions. Lifecycle inventories can include improved esti¬ 
mates non-point source nutrient flows to surface waters by in¬ 
corporating climatic variability. Nutrient discharges to surface 
water are estimated with emission factors as a linear function of 


the annual rainfall rate. Water quality data is required to cali¬ 
brate this model for a given region. In comparison with a stan¬ 
dard approach that uses an average emission factor, the model 
presented here is superior in terms of capturing the variability 
that is correlated to an increase in the size of the hypoxic zone 
in the Gulf of Mexico. 

Recommendations and Perspectives. Lifecycle inventories quan¬ 
tifying nutrient discharges from corn-soybean production should 
include the variability in these flows that occur due to climatic 
conditions. Failure to do so will reduce the LCA's capability of 
quantifying the very significant eutrophication and hypoxia im¬ 
pacts associated with wet years. 

Keywords: Biofuels; corn; eutrophication; hypoxia; non-point 
source pollution; nutrient leaching; soybean; water quality 


Introduction 

Fertilizers used to increase the yield of row crops used for 
food or biofuels can migrate through the environment and 
potentially cause adverse environmental impacts. Nitrogen 
fertilizers have a complex biogeochemical cycle. Through 
their transformations and partitioning among environmen¬ 
tal compartments, they contribute to eutrophication of sur¬ 
face waters at local and regional scales, groundwater degra¬ 
dation, acid rain, and climate change. Phosphate fertilizers 
have a simpler fate in the environment, although leaching of 
soluble and bound phosphorus is an important contributor 
to eutrophication. 

Row crop production in the Midwestern region of the USA 
contributes the highest fluxes of nitrogen and phosphorus 
to the Mississippi River basin [ 1 ] and is considered one of 
the primary contributors to the growing hypoxic zone in 
the Gulf of Mexico. Through a combination of excessive 
nutrient loads and hydrodynamic conditions, a region along 
the coast of Louisiana that is approximately the size of the 
State of Massachusetts is considered ecologically dead most 
summers [2]. Variability in the nutrient loads to the Gulf of 
Mexico is highly dependent on the annual mean river dis¬ 
charge [11, which in turn is a direct consequence of variabil¬ 
ity in the annual rainfall in the Midwest each year. Indeed, 
numerous field studies [3-7] have quantified the substantial 
variation in nutrient leaching during years of light versus 
heavy rainfall. The higher nutrient loads observed during 
wet years are also somewhat correlated to the size of the 
hypoxic zone in the Gulf of Mexico [2,8]. 

Highly sophisticated hydrologic and agricultural models can 
be used to estimate the fate of nitrogen and phosphate nu¬ 
trients in row crop systems and their subsequent fate down¬ 
stream in surface water bodies. Mechanistic models for nu- 
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trient flows attempt to explicitly model the non-linear inter¬ 
relationships between these flows and nutrient use, crop 
productivity, weather patterns, and geographic conditions. 
Application of this type of model requires a vast set of site- 
specific parameters to quantify all of these interrelationships. 
SWAT (soil water assessment tool), for example, is one of the 
most comprehensive watershed models currently available [9], 
especially when linked to a model such as EPIC (erosion pro¬ 
ductivity impact calculator) that estimates field-scale crop 
growth, N and P cycling, runoff and erosion [10]. SWAT has 
been applied to watersheds in central Iowa to predict the 
effectiveness of changes in agricultural practices on the re¬ 
duction in erosion and nitrate loads to surface waters [11]. 
Models such as this require short time steps to delineate 
erosion and nutrient flows caused by individual rainfall 
events because even short term rain events can cause much 
longer-term discharges of nutrients to surface water bodies 
[12-13]. This degree of detail is necessary if the environ¬ 
mental impacts of concern also occur at the same time scale. 
Eutrophication and hypoxia, however, have been correlated 
to seasonal or annual cumulative nutrient discharges [14]. 

Empirical emission factor approaches provide an alterna¬ 
tive means of estimating non-point source nutrient flows as 
a fraction of the fertilizer applied or other measure of an¬ 
thropogenic inputs [15-18]. For example, the International 
Panel on Climate Control (IPCC) [19], which focuses on air 
emissions, recognized the need to quantify leaching to sur¬ 
face waters due to the subsequent denitrification and release 
of N 2 0. The IPCC suggests that nitrogen leached to surface 
waters could be estimated as 30% of the applied fertilizer 
that remains after volatilization, although they recognize 
values could range from 10-80%. The greenhouse gas and 
air emission model (GREET) [15], used leaching rates mea¬ 
sured at several different field test plots [e.g., 20-21] to esti¬ 
mated that 24% of the fertilizer nitrogen (FN) leaches to 
surface water. Neither of these widely used models incorpo¬ 
rates variability in leaching rates as a function of precipita¬ 
tion or other climatic variation. Emission factors are easy to 
apply, and could be appropriate for simple systems such as 
phosphate discharges to surface water. However, they are 
less likely to capture complex interactions among various 
transformations and partitioning of nitrogen among envi¬ 
ronmental compartments, especially as a function of the 
variety of agricultural practices and geologic, hydrologic and 
weather conditions among different farming regions. 

Integrating non-point source nutrient loads associated with 
agricultural activities into lifecycle inventory and impact 
assessments requires a balance between accuracy and rea¬ 
sonable simplicity. Mechanistic models such as those de¬ 
scribed above have not been used in agricultural life cycle 
assessments to date. This is most likely due to the extensive 
parameterization required (soil types, weather, agricultural 
practices etc.). On the other hand, emission factor approaches 
can fail to estimate the range of impacts associated with 
climatic and regional variability. Miller et al. [17] recently 
employed Monte Carlo techniques to provide statistical dis¬ 
tributions of non-point source nitrogen inventory flows. This 
technique includes variability and uncertainty in each indi¬ 
vidual emission factor, but does not integrate the more di¬ 
rect functional relationships among the various nitrogen 


transformation and transport processes. The work presented 
in this paper provides another relatively simple empirical 
approach to integrating variability of nutrient discharges to 
surface water and the resulting eutrophication impacts. 

1 Goals and Scope 

With an increase interest in the use of corn and soybeans for 
fuel production, it is important to understand the relative 
environmental benefits and deleterious impacts associated 
with this growing market. A team of researchers lead by the 
U.S. DOE's National Renewable Energy Laboratory (NREL) 
completed life cycle assessments (LCA) for biodiesel [22] 
and ethanol from com stover [23]. The stover report fo¬ 
cused on the green house gas emission benefits associated 
with biofuels as balanced by potential detriments to soil 
health. Eutrophication was identified as an important issue 
by stakeholders, but limited resources prevented this envi¬ 
ronmental impact category from being addressed. 

The overall goal of this research was to provide a tool to en¬ 
able improved estimates of non-point source nutrient flows from 
row crop production that could be integrated into biodiesel or 
ethanol biofuel lifecycle assessment studies. The objectives of 
the research presented in this paper were to develop, calibrate 
and validate a simple emission factor-based means to accu¬ 
rately quantify nutrient loadings to surface water that reflects 
annual variations due to weather. This methodology was then 
used to quantify surface water quality impacts associated 
with nutrient fertilizers used for energy crop production. A 
set of watersheds that comprise most of eastern Iowa was 
chosen as the geographic region. Ample data describing corn- 
soybean agricultural practices in this region and nutrient 
loadings to the Mississippi River enabled the development 
and calibration of the model for this particular region. 

Applying an agricultural LCA to a specific geographic loca¬ 
tion provides increased accuracy by using site-specific data. A 
balance must be made, however between a site that is so small 
that it is not representative of the broader system versus a 
region that is so large that input data and nutrient flows vary 
substantially over the region and are difficult to characterize. 
The eastern Iowa region selected for this study (Fig. 1) encom¬ 
passes approximately 50 thousand km 2 . It was chosen based 
on its high productivity of corn and soybeans, the availability 
of agricultural practice, yield, and erosion data, and its align- 



Fig. 1: Geographic scope of study area within the State of Iowa includes 
discharge of three major watersheds (Skunk, Iowa and Wapsipinicon) to 
the Mississippi River 
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ment with three major watersheds for which several years 
of water quality data exist. These data were essential for the 
development and calibration of nutrient leaching models. 

County boundaries were used to define the geographic sys¬ 
tem boundary to approximate the overall area encompassing 
the Wapsipinicon, Skunk, Iowa and Cedar River watersheds. 
There is over one order-of-magnitude variation in stream flow 
among the 21 individual years considered (1978-1998), with 
the lowest flows in 1989 and the highest in 1993, correspond¬ 
ing to annual rainfall extremes. Data from 1988-1998 were 
used to develop and calibrate the leaching model, with data 
from 1978-1987 used to validate the overall process. 

The geology, hydrology and land use practices in this region 
are described extensively by Becher et al. [24-25]. Overall, 
approximately 64% of the region is used for row crop pro¬ 
duction. This is mostly dedicated to corn-soybean rotations, 
although some farms, especially in the northeast and east 
central regions, practice corn-corn-soybean rotations [26]. 

Data describing agricultural practices and yields were readily 
available through the USDA National Agricultural Statis¬ 
tics Service (USDA NASS) electronic reports and databases. 
The inclusion of several years in this study allowed variabil¬ 
ity of crop yields and nutrient flows as a function of rainfall 
to be incorporated. The use of average data would have lost 
the important environmental impacts associated with ex¬ 
treme climatic conditions. 


2 Methods 

2.1 Quantification of actual nutrient loads to the 
Mississippi River 

The eastern Iowa watershed system was selected as a geo¬ 
graphic boundary due to the extensive data for total nitro¬ 
gen (TN) and total phosphorus (TP) loads to the Mississippi 
River. The availability of these data provides an excellent 
means of calibrating and validating the nutrient flow model 
against real data. Leaching from fields and in-stream deni¬ 
trification processes were both considered in this calibra¬ 
tion process. 

Data sets quantifying the TN and TP load were available 
from the USGS National Ambient Water Quality Assess¬ 
ment (NAWQA) [25] and the NOAA/USGS Hypoxia in the 
Gulf of Mexico study [27]. Neither of these data sets pro¬ 
vides both TN and N0 3 -N loads for all three watersheds of 
interest over the entire time period, so some of the loads 
were estimated based on statistical correlations. TN and TP 
loads are presented in Table 1. Values for the Cedar River - 
which is included within the overall Iowa River watershed - 
are included because they provided the best correlation for 
the fraction of TN load associated with the Wapsipinicon 
River. For the three years that a complete data set for all 
four watersheds was available [25], it was determined that 
the Wapsipinicon TN loads were approximately one-third 


Table 1: Actual nutrient loads to the Mississippi River from E. Iowa watersheds 


Year 

Total Nitrogen (t/y) 

Total Phosphorus (t/y) 

Wapsipinicon 

Iowa 

Skunk 

Cedar 

Wapsipinicon 

Iowa 

Skunk 

Cedar 


Data used to validate model 


1978 

6,990 

58,332 

26,469 

19,970 

319 

2,661 

1,751 

700 

1979 

18,031 

92,205 

29,819 

51,516 

490 

4,252 

1,910 

1,679 

1980 

11,185 

45,342 

11,098 

31,957 

226 

1,981 

638 

992 

1981 

7,969 

42,255 

11,884 

22,767 

201 

1,836 

668 

718 

1982 

17,255 

115,248 

41,275 

49,300 

523 

4,509 

2,381 

1,478 

1983 

24,964 

120,828 

38,033 

71,325 

542 

4,641 

1801 

2,241 

1984 

17,737 

95,607 

35,955 

50,678 

459 

3,792 

1,925 

1,623 

1985 

6,290 

44,659 

19,346 

17,970 

191 

1,702 

827 

522 

1986 

15,672 

103,254 

33,968 

44,778 

443 

3,811 

1,961 

1,324 

1987 

3,865 

36,485 

1,6045 

11,044 

154 

1,464 

638 

360 


Data used to calibrate the model 


1988 

2,200 

20,179 

5,087 

6,290 

81 

913 

146 

236 

1989 

1,227 

7,616 

2,189 

3,508 

52 

557 

118 

155 

1990 

9,514 

69,219 

22,483 

27,196 

354 

3,108 

1711 

844 

1991 

19,618 

96,879 

22,729 

56,079 

416 

3,777 

1112 

1,778 

1992 

13,841 

88,422 

23,860 

39,564 

326 

2,893 

1,295 

1,027 

1993 

35,784 

199,911 

51,287 

102,290 

1,158 

10,341 

4,673 

3,520 

1994 

9,336 

46,969 

8,026 

26,688 

172 

1,711 

365 

683 

1995 

8,284 

54,321 

15,877 

23,680 

248 

2,177 

1,140 

653 

1996 

9,990 

57,600 

30,000 

31,900 

540 

2,360 

3,960 

1,560 

1997 

21,700 

75,100 

23,500 

52,200 

440 

3,120 

990 

1,860 

1998 

31,100 

154,000 

49,400 

97,000 

476 

3,930 

4,140 

2,470 


1996-1998 data from Becher et al. [24] 

1978-1995 data for Skunk, Iowa and Cedar Rivers from NOAA/USGS Hypoxia study 
1978-1995, Wapsipinicon R. - Estimated based on 0.35 * TN ( cedar), 0.062 * TP (C edar+skunk+iowa) 
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1960 1970 1980 1990 2000 

Fig. 2: Historical nutrient use of commercial nitrogen (FN) and phosphate 
(FP) fertilizers for corn and soy in Iowa: a) application rate; and, b) fraction 
of fields applying fertilizer 

(0.35T0.06) of the Cedar River values. This fraction was 
used to estimate the total nitrogen load to the Mississippi 
River during the period 1978-1995. TP values for the 
Wapsipinicon (1978-1995) were estimated as a fraction 
(0.062T0.015) of the total loads from other rivers. 

Becher et al. [24] provide an estimate of nitrogen sources in 
Eastern Iowa. They show that a combination of fertilizer 
nitrogen (FN) and animal manures contribute over 90% of 
the TN discharged to the Iowa and Skunk Rivers. Thus, we 
assumed here that the inclusion of animal manures as a TN 
source would enable us to perform the most accurate cali¬ 
bration. FN used in this region was calculated from pub¬ 
lished county totals for corn and soy acreage, statewide av¬ 
erage FN application rates, and the fraction of farms applying 
FN to their crops (Fig. 2). Animal manure generated within 
the region was estimated based on county level data for live¬ 
stock inventory, which were available from the agricultural 
census on-line for the years of 1988,1989,1990,1992,1997 
and 2002 [28]. Cattle, cows, hogs and sheep were included 
in the manure calculations. Hogs were observed as having 
the greatest contribution of manure nitrogen. County aver¬ 
ages values were used for the missing years and all years of 
the validation period (1978-1987). Nutrient composition 
(N,P) for each type of animal manure was obtained from 
the USDA Agricultural Waste Field Handbook [29]. 


the erosion potential [30]. Finear, exponential and other non¬ 
linear equations have resulted from data fitting exercises. In 
all cases, annual rates or flows were deemed most suitable 
for these empirical relationships. 

Bakhsh et al. [7] found that rainfall was the most significant 
variable affecting nitrate leaching from a research field in 
eastern Iowa. They could describe 79% of the total vari¬ 
ability in nutrient loads to surface water as a linear function 
of the annual subsurface drainage volume. Drainage, in turn 
was linearly related to rainfall. 

Review of the TN and TP loads from eastern Iowa to the 
Mississippi River can help to identify which of these poten¬ 
tial variables could be used to provide accurate nutrient load 
predictions while maintaining a reasonable level of simplic¬ 
ity. Surprisingly, the TN and TP loads from eastern Iowa 
watersheds to the Mississippi River were highly correlated 
(1978-1995 data for Skunk, Iowa and Cedar Rivers from 
NOAA/USGS Hypoxia study, Table 1). A linear regression 
shows that TN loads for these three rivers are 23.3 times 
higher than TP loads (R 2 =0.88). This indicates that, even 
with the greater solubility and much more complex trans¬ 
formation processes possible with nitrogen species, it be¬ 
haves in a similar fashion to phosphate species. The TN and 
TP loads for these three rivers are also highly correlated to 
the average annual river flow rate (R 2 =0.90, 0.91, respec¬ 
tively). Given the preponderance of evidence that water trans¬ 
port processes dominate the delivery of fertilizer nutrients 
to surface water bodies, it was hypothesized that a simple 
linear function that predicts nutrient loading as a function 
of rainfall - used here as a surrogate for river flowrates - 
would capture the predominant mechanisms controlling 
these average annual flows and enable accurate prediction 
of TN and TP loads to the Mississippi River. 


2.3 Mathematical form for the nutrient leaching model 

Nutrient flows to the Mississippi River were calculated with 
a standard emission factor approach with the adaptation 
that the emission factors varied as a function of rainfall. 
The general leaching models for TN and TP are defined based 
on a fraction of the applied nutrient load. The nitrogen model 
is amended to account for the subsequent loss of nitrate via 
denitrification in tile drains and local streams. 


t n - 
l MS - 


I 

j=county 


e+L 


N 




+ L 


N 





TN,N03 


' fde. 


( 1 ) 


2.2 Basis for nutrient leaching model 

Previous research efforts to establish empirical formulas to 
estimate nutrient loads or fluxes to the Mississippi River 
have incorporated a wide variety of parameters. Some ex¬ 
amples include fertilizer application rates in the current and 
preceding years [16,30,31], total anthropogenic nitrogen 
inputs [32], residual nitrogen (inputs minus outputs) during 
the previous year [31], average river discharge rates or stream 
yield (flow/ watershed area) [31-33], and for phosphorus, 


rP 

l ms 


I 

j=county 



+ L 



( 2 ) 


Where F is the mass load (t y 1 ) for the region with super¬ 
scripts defining the nutrient (N or P) and subscripts defining 
the mass discharged to the Mississippi (MS) or added to the 
system as commercial fertilizer (f) or animal manure (an). 
The f terms are emission factors that define fractions of the 
total flows. Subscripts define the fraction of the applied nu- 
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trient leached to the surface water (SW), the fraction of the 
total nitrogen that is nitrate (TN,N03), and the fraction of 
nitrate in the surface water that is denitrified within the 
smaller rivers before discharge to the Mississippi River (de,ri). 
Terms with subscript j are quantified at the county level. 

The commercial fertilizer loads for each county (L f k ) were 
estimated based on appropriate historical values for the area 
planted (A ; ) (ha) multiplied by the fertilizer application rates 
(N fi ) (t ha -1 y 1 ) and fraction of hectares fertilized (f fi ). 


Pj ’^f,c' k f)c + A soy,pj '^f,soy’%,soy 


( 3 ) 


L f. A c, 

Subscripts c and soy refer to corn and soybeans, respectively 
and superscript k refers to nitrogen or phosphorus. Data for 
the area planted in each crop were available for each year at 
the county level. Fertilizer application rates were available 
only at the state level and were assumed to be representative 
for the eastern region. These data, acquired from standard 
agricultural databases, [34-35] are summarized in Fig. 2. 


Eq. (1) and (2) also require estimates of the fraction of ni¬ 
trogen or phosphorus fertilizer that leaches to the surface 
water for each county (f sw k ). These terms were estimated as 
a linear function of rainfall: 


f s k w =b k +m k -Rj (4) 

Where m and b are the slope and intercept of a linear func¬ 
tion, and Rj is the annual rainfall (mm), determined at the 
centroid of each county [36]. The linear coefficients in Eq. (4) 
and denitrification rates (f deri ) (Eq. (1)) were calibrated by 
minimizing the difference in the cumulative mass loading 
between the nutrient leaching model predictions (Eqs. (1)- 
(4)) and actual TN and TP loads to the Mississippi over the 
11 year period used for to calibrate the model (1988-1998). 
Two key assumptions were made in this analysis: 

1. The same fraction of the nutrients leaches from animal 
manure and inorganic commercial fertilizer. This assump¬ 
tion is also use in the IPCC nitrogen flow model. 

2. Denitrification of N0 3 in rivers also varies some with 
rainfall [37]. The fraction of the TN that is N0 3 was 
required to estimate denitrification. This fraction was 
defined as an average for watershed discharge data avail¬ 
able 1988-1998, with a weighted average between the 
three rivers (f TNN03 = 0.83). 

2.4 Allocation of nutrient leaching flows 

In an LCI, it is important to not only quantify flows through 
the environmental compartments, but also to allocate specific 
fractions of the total flows among the various products and 
by-products generated. This is particularly difficult in an agri¬ 
cultural system when there are synergistic affects associated 
with crop rotation. Van Zeijts et al. [38] address this issue for 
allocating fertilizer among various crops in the Netherlands. 
Their overall approach is used here as well. They assume: 

• All nitrogen can be allocated to the crop for which it is 
applied in the year that it is applied 

• All P and K fertilizers should be allocated among the crops 
over the entire crop rotation in proportion to the fraction 
of nutrient that is harvested with each of the grains. 


They justify the differences in these approaches based on 
the intent of the fertilizer application and its environmental 
fate. Nitrogen, for example, is typically applied separately 
for each crop in proportion to that crops needs. Excess fertili¬ 
zation is avoided due to its likely loss to leaching. Phosphate, 
however, behaves differently. Excess nutrients applied as fer¬ 
tilizer for one crop is conserved in the soil over a longer pe¬ 
riod. It is often applied only to one crop with the intent that it 
be consumed by the plants over the entire rotation. The con¬ 
sistency in the percentage of fields onto which nitrogen fertil¬ 
izer (FN) and phosphate fertilizer (FP) are applied (see Fig. 2) 
illustrates that indeed, many farmers chose not to apply FP to 
fields in a soy rotation year, even though the soy beans need 
phosphorus for growth. These farmers rely on FP remaining 
in the soil from the previous year. 

The suitability of the Van Zeijts et al. [38] approach for 
nitrogen allocation between corn and soybeans could be 
argued. Nitrogen fixation by soy increases the nitrogen con¬ 
tent in the soil and the soil mineralization rate during soy 
and corn years, thus providing extra nitrogen that would not 
otherwise be available [39]. Likewise, excess nitrogen applied 
to the soil during a corn year and immobilized with in soil 
organic matter could become available to soy plants in the 
following year [40]. The 'carry over' of nitrogen fertilizer is 
indeed observed in leaching rates that are generally no dif¬ 
ferent on a kg/ha basis in years that corn is grown versus 
years that soy is grown [5,7]. Thus, the approach used by 
van Zeijts et al. [38] for nitrogen allocation is simplistic. 

An alternative approach for allocating the amount of fertil¬ 
izer nitrogen (FN) leached from corn and soy fields can be 
based on field data that quantifies actual leaching rates from 
corn fields versus soy fields. Data from each of the studies 
at the Nashua IA agricultural testing facility [5,7,41] indi¬ 
cate that even though little to no FN is applied to soy fields, 
nitrate fluxes (t ha -1 y 1 ) in a given year are approximately 
equal from each type of field. Statistical analysis showed 
that there are no trends in the fraction of nitrate that leaches 
from corn fields as a function of year or tillage practices 
(a=0.05). Thus, all data were lumped together (n=42) to show 
that 51.1% (standard deviation = ±7.2) of the total amount 
of nitrate that leached on a per hectare basis could be allo¬ 
cated to the corn. The total mass of nitrogen discharged to 
the Mississippi River from commercial fertilizers used in 
eastern Iowa that can be allocated to corn (L N MSc ; t y _1 ) can 
be determined as an area weighted average, with slightly 
higher losses from corn fields: 


t N 

l MS,c 


2 A c, Pi 
j=county 


x < 

j=county 


i 

j=county 


A c,P, + 


1 fN -leach,, 


■A 


s °y<Pi 


v /n 

L r ' f SW. ' “ fTN,N03 ' fde,ri ) 


( 5 ) 


where the fraction allocated to corn (f N _i eac h, c ) was set at 51%. 
The summation on the right represents the commercial fer¬ 
tilizer component of Eq. (1). 
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3 Results and Discussion 
3.1 Nutrient leaching model 

The nutrient leaching models for TN and TP (Eqs. (1)—(4)) 
were calibrated independently. The calibrated model esti¬ 
mates are within 1.1% of the actual cumulative TN load 
and 0.3% for the total mass of TP discharged over the 11 
years. Values for the calibrated model parameters used are 
included in Table 2. The model estimates for L N MS and L P MS 
were very sensitive to the slope, intercept and rain data used 
in the regression equations. Very small modifications in the 
slope and intercept helped substantially to improve the qual¬ 
ity of the fit. Application of the high denitrification rates 
(f d e ri) suggested by deVries et al. [16] for the Netherlands 
resulted in very low estimates of nutrients discharged to the 
Mississippi (L MS ). Thus, lower denitrification rates were re¬ 
quired in both high and low flow years to calibrate the model. 
The TN model was not very sensitive to the fraction of the 
total nitrogen that is in the form of nitrate (f TN nos), so this 
parameter was not adjusted. 


Table 2: Calibrated emission factors used in the nutrient leaching models 


Model 

Calibrated 

Comments 

Parameter 

Value 



Nitrogen leaching model (Eq. (1), (4)) 


Intercept 

-0.545 

Highly sensitive to small changes 

Slope 

0.000975 

Highly sensitive to small changes 

fTN,N03 

0.83 

Not sensitive 

fde.ri 

0.2 

For years with high flow - highly sensitive 

fde.ri 

0.3 

For years with low flow. (<900 mm) 

Phosphorus leaching model (Eq. (2) and (4)) 

Intercept 

-0.11 

Highly sensitive to small changes 

Slope 

0.000174 

Highly sensitive to small changes 


Fig. 3 presents the model predictions and actual TN and TP 
loads to the Mississippi River. The low TN and TP loadings 
in 1988 and 1989 correspond to drought years and the very 
high loadings in 1993 correlated to heavy rains and flood¬ 
ing. With a few exceptions, the calibrated model represents 
both the trends and the absolute value of the TN and TP 
loads. For individual years, the greatest percent errors in 
TN loads were for 1988 (-88%) and 1989 (+160%). 1988 
resulted in the greatest error for TP loads as well (-99%), 
although it is apparent that the absolute errors in these years 
were small. Calibrated estimates of TN loads for seven of 
the years were within 25% of the loads determined from 
water quality data. 

The calibrated model was applied to the 1978-1987 period 
for validation (see Fig. 3; open circles). Both the TN and TP 
predicted loads were very representative of the trends and value 
of the actual loads to the Mississippi River. For TP loads, the 
absolute value of errors in the annual loads ranged from 0.15% 
(1983) to 104% (1985), with the error for the cumulative 11 
year load 20%. Similar results were obtained for the TN model. 
The maximum error in the annual load was 95% (1981). 
Both the fact that predictions for seven of the ten years were 
within 25 % of the actual data and the ten-year cumulative 
load within 14% of the actual load verify the validity of this 




Fig. 3: Actual and predicted loads of total nitrogen (a) and total phospho¬ 
rus (b) from eastern Iowa watersheds to the Mississippi River. Solid circles 
represent the calibrated model and open circles show the model validated 
against data not used for the calibration process 

relatively simple approach for estimating nutrient loading 
to surface water from fertilizer use on row crops. 

The quality of the nitrogen leaching model presented here 
can also be assessed by comparing the emission factor for 
leaching (P^y) and the nitrogen yields to those measured by 
others. Over the entire region, the model presented here es¬ 
timates the load of TN to surface waters through tile and 
base flow to be 6.4-85% of the applied fertilizer. Distribu¬ 
tions of these emission factors at the county level for both 
nitrogen and phosphorus are presented in Fig. 4, along with 



0.0 0.2 0.4 0.6 0.8 1.0 

Fraction Applied Fertilizer to SW 

Fig. 4: Histograms illustrate the distribution of county-level rainfall (a) and 
the fraction of nitrogen and phosphate fertilizers that are delivered to sur¬ 
face water bodies 
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Table 3: Range of estimates - emission factors for quantifying fraction of 
commercial nitrogen fertilizer that leaches with water 


Emission 

Factor 

Comments, references 

0.316 

Median - this study (Fig. 4) 

0.24 

Used in LCA for subsequent denitrification and N 2 0 
emissions for transportation fuels (GREET [15]) 

0.30 

Global average value used for subsequent 
denitrification and N 2 0 emission (IPCC [19]) 

0.25 

(0.015-0.64) 

Measured at downstream location. Estimated from 
cumulative load in Raccoon River (central IA) over 

11 years versus cumulative FN used. [3] 

0.36 

(0.08-0.59) 

Measured - 3 year field scale study. Nashua IA [4] 

0.26 

(0.05-0.76) 

Measured - 6 year field scale study. Nashua IA 
(NE IA). Rainfall rates most significant variable [5] 

0.42 

(0.09-0.85) 

Measured - Walnut Creek watershed study, central 
IA, average of several drainage-scale systems [ 6 ] 

0.47 

Overall 6 year cumulative average fraction FN 
leached from single C-S rotation field in Walnut 
Creek watershed, central IA (1992-1997) [ 6 ] 


the correlated distribution of rainfall. As expected based on 
the lower solubility of phosphates, a much smaller fraction 
of phosphate fertilizer enters the surface water than of ni¬ 
trogen fertilizer. The median of emission factors for nitro¬ 
gen (0.32; Table 3) is within the range of averages presented 
by others (0.24-0.47). The skewed distribution for nitrogen 
leaching shown in Fig. 4 exemplifies the importance of pre¬ 
dicting much higher nutrient loads during rainy years. 

When the amount of nitrogen leached is presented as a yield 
(kg ha -1 y -1 )> predicted TN loads are also quite comparable 
to other studies from the Midwestern United States (Table 4). 
Given the consistency of the yields predicted with the nitro¬ 
gen model used here and those measured by others, and the 
accuracy of the validated TN leaching model (see Fig. 3), 
the estimated discharges of TN from row crops and to the 
Mississippi River are considered realistic. 

The model estimates of the fraction of the fertilizer phos¬ 
phorus that leaches, -0-12% (median 4.2%, see Fig. 4), 
compares quite well with global estimates of phosphorus 
fertilizer leaching rates (1-5%) [42] and those measured in 
a limited number of field studies in Iowa (5-15%) [43]. 


Table 4: Comparison of nitrogen yields estimated with the leaching model 
and those measured by others 


Range of Values 
(kg-N/ha crop/yr) 

Comments, references 

5.8-64 

TN - This study (median: 32 kg/ha) 

38-64 

NO 3 -N in tile drain water, IL 30-ha field [51] 

11-107 

NO 3 -N in tile waters, 4 different tillage systems 
in IA[41] 

4-46 

NO 3 -N in tile waters, field study in IA [7] 

2-60 

NO 3 -N in tile waters, 8 year watershed study, IA [33] 

10-80 

TN from OH watershed [52] 

38±16 

Model results for corn in Midwest USA [17] 

22±13 

Model results for soybeans in Midwest USA [17] 


The approach used by LCI models such as GREET [15] (TN 
leaching = 0.24 • L f ) provides an excellent estimate of the TN 
load for 'average' years (see Fig. 3). However, the use of a 
constant emission factor does not capture the wide variability 
in leaching rates that are correlated to rainfall. Errors with the 
GREET approach are as high as 900% (1989), with an over¬ 
all average and standard deviation of errors over the 21 year 
period of -66%±211%. In contrast, the average error with 
the variable leaching model developed here is -17%±53%. 
The biggest concern with the use of a constant emission factor 
for nitrogen leaching is the inability of this approach to pre¬ 
dict the extremely high nutrient loadings that occur in rainy 
years. The actual loads in 1993 were over three times the value 
predicted using the GREET emission factor. The extreme nu¬ 
trient loads that occurred in 1993 led to a record size of the 
hypoxic zone in the Gulf of Mexico (18,000 km 2 ) [2]. The 
ability to accurately predict extreme LCI nutrient loads with 
the model proposed here will also lead to an improved abil¬ 
ity to predict their related environmental impacts. 

3.3 Allocating nutrient flows between corn and soybeans 

The allocation process is an important step to define of how 
much of the total environmental emissions can be attributed to 
corn versus soybean. Estimates of the TN leached that can be 
allocated to corn between the two methods considered here 
ranged from 60±4% (Eq. (5)) to 98±1% (Van Zeijts et al. [38] 
method). The balance of the TN load is allocated to soybeans. 
The difference between these methods stems from a poorly un¬ 
derstood symbiosis of nitrogen flows within the corn-soybean 
rotation that is difficult to integrate into an LCA. The alloca¬ 
tion of phosphorus between crops shows that 61±4% should 
be allocated to corn in this system. The standard deviations 
in the above numbers represent variability between different 
years. Unlike the total load of nutrients discharged, there is 
little variability in the fraction of nutrients allocated between 
the crops among years with variable rainfall. 

3.4 Eutrophication impact 

In life cycle assessments, equivalency factors are often used to 
aggregate nitrogen and phosphorus species to quantify the overall 
eutrophication potential. The Redfield molar ratio (16:1) for 
the N:P content in aquatic cell mass is used here to quantify 
the equivalency factors [44]. In reality, water bodies have a 
broad range of ratios in the limiting N:P ratios that enable 
excess algal growth and eutrophication [45]. For example, a 
higher N:P ratio of 22:1 can apply to eutrophic waters [46]. 

The eutrophication potential (t y _1 as N0 3 -N equivalents) 
discharged to local streams and then the Mississippi River 
includes total nitrogen and total phosphorus loads predicted 
from the calibrated leaching models. This potential is com¬ 
pared to several different benchmarks to define the poten¬ 
tial impact of these discharges. Bench marks include: 

• Median background concentration of TN (0.13 mg/L) 
and TP (0.021 mg/L) in streams and rivers in the corn 
belt region and great plains [47] 

• Proposed water quality standards for TN (2.18 mg/L) 
and TP (0.0763 mg/L) in the corn belt region [48,49] 

• Recommended 30% reduced load of 1980-1996 average 
TN load to the Gulf of Mexico to reduce hypoxia [50] 
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Fig. 5: Predicted eutrophication impact based on calibrated leaching mod¬ 
els compared with various benchmarks 


Concentration limits were multiplied by the actual annual 
mean stream flow rate and converted to an eutrophication 
potential (TN+TP, as N0 3 -N) to determine the load. 

The estimated eutrophication potential for the period 1988- 
2000 is compared to the applicable benchmarks in Fig. 5. 
The total modeled eutrophication potential is much higher 
in all years than the load estimated from proposed water 
quality standards. This indicates that unacceptably high lev¬ 
els of biomass growth, decay, and oxygen depletion will 
occur in the streams and rivers in eastern Iowa. Utilizing 
the higher N:P ratio of 22:1 for eutrophic waters would 
show greater contributions from TP and even higher total 
eutrophication potentials. 

As the nutrients are transported downstream in the Missis¬ 
sippi River, they have additional eutrophication consequences 
in the Gulf of Mexico where they contribute to hypoxia. 
The goal of a 30% reduction in the average TN load estab¬ 
lished by a federally endorsed task force on hypoxia [50] is 
shown as a straight line at -500,000 t y 1 as N0 3 -N. Mod¬ 
eling work by Scavia et al. [14] suggest that the goal set by 
Brezonik et al. [50] is too low and that a 40-45% reduction 
in average TN loads would be required to limit the size of 
the hypoxic zone. The 45% reduction in average TN loads 
would lower the acceptable eutrophication potential to 
421,000 t y 1 as N0 3 -N. 

The eastern Iowa contributions to the total eutrophication 
potential exceed the 30% reduction goal in 5 of the 13 years 
shown in Fig. 5. In contrast, the use of a constant emission 
factor for fertilizer leaching (24%) to estimate the nitrate 
load to surface water would predict that the eutrophication 
potential discharged to the Mississippi River is always less 
than the suggested -500,000 t y 1 N0 3 -N target. This ex¬ 
ample clearly shows the necessity of incorporating the cli¬ 
matic variability in loads on actual environmental impact. 

4 Conclusions, Recommendations and Perspective 

The nutrient flow model was calibrated and verified for a 
set of watersheds covering most of the eastern half of Iowa 
for a 21-year period that includes both drought and flood 
years. Water quality data were available for these water¬ 


sheds allowing the TN and TP loads to be calibrated to ac¬ 
tual measured discharges. The incorporation of annual vari¬ 
ability in nutrient loads due to rainfall has not been used 
before with LCAs for non-point source pollutants and pro¬ 
vides an excellent approach for quantifying the true vari¬ 
ability. Although this approach helps to improve the esti¬ 
mates and variability in non-point emissions from nutrient 
leaching, it also makes the model and approach very site 
specific. Thus, the specific linear model presented here to 
estimate the fraction of fertilizer that leaches should not be 
applied directly to other locations. For example, western 
Iowa has sandier soil, which would allow a greater mass of 
nutrients to infiltrate to groundwater, with less to surface 
water. Although it is not intended that the site specific val¬ 
ues presented here be used elsewhere, the general approach 
and framework is valuable and can be applied to other sites 
that have suitable water quality data. 

The results of this analysis show that the eutrophication 
potential for the current system already exceeds acceptable 
limits. TN and TP discharges from row crops in eastern IA 
exceed the maximum loads defined by the proposed water 
quality standards every year. Limits established by the goal 
of a 30% reduction in the average TN load to reduce hy¬ 
poxia are also exceeded in approximately half of the years. 

Lifecycle inventories quantifying non-point source nutrient 
discharges from corn-soybean production should include the 
variability in these flows that occur due to climatic condi¬ 
tions. Failure to do so will reduce the LCA's capability of 
quantifying the very significant eutrophication and hypoxia 
impacts associated with wet years. 
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